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It is well known that the kinetic exchange interaction between single-occupied magnetic orbitals 
(s-s) is always antiferromagnetic, while between single- and double-occupied orbitals (s-d) is always 
ferromagnetic and much weaker. Here we show that the exchange interaction between strongly 
anisotropic doublets of lanthanides, actinides and transition metal ions with unquenched orbital 
momentum contains a new s-d kinetic contribution equal in strength with the s-s one. In non- 
collinear magnetic systems, this s-d kinetic mechanism can cause an overall ferromagnetic exchange 
interaction which can become very strong for transition metal ions. These findings are fully con¬ 
firmed by DFT based analysis of exchange interaction in several Ln 3+ complexes. 


I. INTRODUCTION 


Anderson’s kinetic exchange interaction [1, 2] is ubiq¬ 
uitous in magnetic molecules [3, 4] and insulating ma¬ 
terials [5, 6]. In particular, the kinetic mechanism has 
been found as dominant contribution to the exchange in¬ 
teraction in various transition metal compounds. The 
mechanism has been also often advocated as reason for 
orbital ordering in transition metal oxides with orbitally 
degenerate metal sites [7], especially, in magnetoresistive 
manganese oxides [8]. 

In all these cases the magnetic orbitals are real and the 
exchange interaction in the case of non-degenerate sites 
is described by Heisenberg Hamiltonian, H = — j • S 2 . 
The kinetic exchange interaction originating from virtual 
electron transfer between single-occupied orbitals (s-s) is 
always antiferromagnetic [1]: 


Js-s 




(i) 


where t is the transfer parameter and Uij is the electron 
promotion energy from site i to site j. 

On the contrary, the electron delocalization between 
double-occupied and single-occupied orbitals (s-d) al¬ 
ways results in a ferromagnetic contribution (the Good- 
enough’s mechanism [6]): 
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where Jr is the Hund’s rule coupling constant and t is the 
transfer integral between corresponding orbitals. Given 
the typical ratio Jr/U ~ 0.1 [1, 2], J s _d is by one order 
of magnitude smaller than J s _ s . Then for comparable 
electron transfer parameters in s-s and s-d processes, the 
overall coupling is antiferromagnetic, J = J s -s~\-Js-d < 0. 

A weak ferromagnetic interaction (2) is observed when 
the electron transfer between the single-occupied orbitals 
is negligible or zero, which is achieved for certain geome¬ 
tries of the exchange bridge [2, 5, 6]. Similar ferromag¬ 
netic contribution appears also for electron delocaliza¬ 
tion between single-occupied and empty orbitals, as well 


as in the case of degenerate magnetic orbitals (Kugel- 
Khomskii model) [7]. In all these cases the ferromagnetic 
kinetic contribution arises in the third order of perturba¬ 
tion theory after t and Jr. 


The kinetic exchange mechanism is equally important 
in / electron systems such as lanthanide and actinide 
compounds [9-12]. However its realization in these ma¬ 
terials is expected to be different from transition metal 
compounds due to a more complex structure of multielec- 
tronic states on the metal sites, involving complex mag¬ 
netic orbitals. The last are stabilized by strong spin-orbit 
coupling in lanthanides and actinides giving rise to un¬ 
quenched orbital momentum in their low-lying multiplets 
[13], which persists in any geometry of their environment. 
Unquenched orbital momentum also occurs in many tran¬ 
sition metal complexes and fragments when the latter 
possess cubic [14] or axial [15] symmetry, and it was 
proved that its effects can persist also under significant 
deformations of the ligand environment [16]. Nonethe¬ 
less, despite numerous examples of strongly anisotropic 
magnetic materials with unquenched orbital momentum 
on the metal sites, the basic features of kinetic exchange 
interactions in them have not been yet elucidated. 


In this work the kinetic exchange interaction for metal 
sites with unquenched orbital momentum is investigated. 
On the example of strongly axial doublet states, we show 
that the paradigm of active magnetic orbitals as always 
belonging to half-filled ones does not hold for strongly 
anisotropic systems with unquenched orbital momentum 
on sites. In such systems the kinetic exchange interaction 
between single- and double-occupied orbitals is found to 
be of equal strength with conventional kinetic exchange 
interaction between single-occupied orbitals and can even 
make the entire interaction ferromagnetic. Contrary to 
the Goodenough’s mechanism (2), the s-d kinetic contri¬ 
bution found here appears already in the second order of 
perturbation theory being of the form (1). 
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FIG. 1. Electron transfer processes between doublets with unquenched orbital momentum, (a) Scheme of / orbital levels in 
axial crystal field. Numbers in the left side stand for the orbital angular momentum projection on the axis of the field. Pictures 
in the right side show the corresponding real orbitals (for m ^ 0, only one of the two is shown). The electron configurations 
correspond to the wave functions of the Kramers doublet for Dy 3+ with maximal projection of J. (b) Electron transfer processes 
between two collinear Dy 3+ ions in axial Kramers doublets with maximal total momentum projection, Mj = ±15/2. s-s and 
s-d processes are shown by solid and dashed lines, respectively, m, n stand for orbital momentum projections on the direction 
of anisotropy axis on each metal site. Right plots correspond to reversed spin configuration on the site 2. 


II. RESULTS 

A. Doublets with unquenched orbital momentum 

Metal ions are often characterized by non-zero or¬ 
bital momentum L [14, 17]. However, in order to keep 
(part of) it unquenched in complexes and crystals, the 
metal ions should also possess strong spin-orbit coupling 
which splits strongly the atomic (ionic) LS term in multi- 
plets corresponding to definite total angular momentum 
J = L±S [17]. This is a standard situation in lanthanides 
and actinides [11, 13]. Transition metal complexes in a 
threefold degenerate orbital state possess an unquenched 
orbital momentum corresponding to an effective L = 1 
[14]. In this case the spin-orbit coupling leads to the for¬ 
mation of multiplets corresponding to total pseudo mo¬ 
mentum j = S ± 1, S , |S — 1|. 

In low-symmetry crystal field, the (pseudo) J- 
multiplets on metal sites split into Kramers doublets in 
the case of odd number of electrons, or into singlets in 
the case of even number of electrons. The singlets in the 
latter case form quasi doublets for large J or perfectly 
degenerate (Ising) doublets in environments of axial sym¬ 
metry [18]. In all these doublets the two wave functions 
are related by time inversion [17]. Besides, they are mag¬ 
netic and contain a significant contribution of orbital mo¬ 
mentum [19]. The latter necessarily implies that the mag¬ 
netic orbitals and the wave functions of the doublets are 
complex. In the following we consider the simplest case 
of an axial crystal field, in which the atomic orbital wave 
functions preserve the projection of orbital momentum 
(1) on the symmetry axis (m). The crystal-field orbitals 
are twofold degenerate with respect to the sign of the pro¬ 


jection m and are described by the eigenfunctions |Z, ±ra) 
(Fig. la). For more than half-filled atomic orbital shell 
l N , N > 21 ± 1, the ground atomic multiplet corresponds 
to J = L ± S while the wave functions corresponding 
to the maximal projection, Mj = ±J, are represented 
by single Slater determinants. An example is the ground 
Kramers doublet of Dy 3+ ion in strong axial crystal field 
shown in Fig. la [18]. Further we consider this kind 
of axial magnetic doublets only [20], which allows us to 
describe the exchange mechanism in the simplest way, 
though the discussed effects are general for all doublets 
with unquenched orbital momentum. It is worth men¬ 
tioning that the doublet states | J, ±J) appear quite often 
in the ground state of lanthanides and represent a great 
interest for the design of single-molecule magnets [21]. 

B. Exchange interaction for collinear doublets 

The kinetic exchange interaction between doublet 
states is conveniently described by pseudospin formal¬ 
ism [17], in which the doublet eigenfunctions | J, ±J) are 
put in correspondence to eigenfunctions 11/2, ±1/2) of an 
effective S = 1/2. First, we consider the case of collinear 
doublets, when their main magnetic axes are parallel. 
Since one-electron transfer processes neither can switch 
nor mix the two doublet wave functions on each metal 
site, for relatively large J, the exchange Hamiltonian re¬ 
duces to the following Ising form [22]: 

H = -JS 1z S 2z , (3) 

where Si Z is the z component of the S^, directed along 
the main magnetic axis on the corresponding metal site 
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[23]. In this case the exchange parameter J is simply 
derived from the difference between energies of antifer¬ 
romagnetic and ferromagnetic configurations in Fig. lb, 
J = 2(Eaf — E f ). We calculated separately the contri¬ 
butions from s-s and s-d processes (Fig. lb) to Eaf and 
E f in the second order of perturbation theory after elec¬ 
tron transfer. This yields the following contributions of 
s-s and s-d processes to the exchange coupling constant 

J : 


JS - S ~ (c/l2 + [/2i) ^ ? |tm ’" n|2 ’ 
Xd = -T- £ Y, (l* m ,-nr-|t m ,n| 2 ) 
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where m and n denote the orbitals on site 1 and 2, 
respectively, by corresponding angular momentum pro¬ 
jections (Fig. lb), and Si and di indicate the sets of 
single- and double-occupied orbitals in the electron con¬ 
figuration | J, J) of site i, respectively. For example, for 
Dy 3+ ion (site 1 in Fig. lb) Si = {—3, —2, —1,0,1} and 
di = {2,3}. In these equations, £ m?n are electron trans¬ 
fer parameters between orbitals m and n. Note that we 
do not include effects oc Jh (Goodenough’s mechanism) 
as being much weaker compared to the s-d contribution 
found here (vide infra). 

While Eq. (4) looks as a standard expression for the 
s-s kinetic exchange parameter [1, 2], the s-d kinetic con¬ 
tribution, Eq. (5), does not appear for isotropic magnetic 
systems in this lowest order of the perturbation theory. 
We can see that it contains electron transfer terms of both 
signs, be., antiferromagnetic and ferromagnetic contribu¬ 
tions. The terms with m = 0 in the two brackets of Eq. 
(5) mutually cancel because of the relation |to, n | = |to,-n| 
[24]. Another evident cancellation occurs for terms with 
n = 0. The other pairs of terms in Eq. (5), with m,n/ 0, 
will not cancel each other unless the metal-ligand-metal 
fragment possesses special point symmetry. Therefore, 
for general geometry of exchange-coupled pairs, the s- 
d kinetic exchange is operative and represents a new 
mechanism of exchange interaction, proper to strongly 
anisotropic metal ions with unquenched orbital momen¬ 
tum only. The peculiarity of this mechanism is that it 
is of the order ~ t 2 /U, z.e., of similar strength as the 
s-s kinetic exchange, Eq. (4). However, at variance with 
the s-s kinetic exchange, the s-d exchange can be both 
antiferromagnetic and ferromagnetic as Eq. (5) shows. 

Due to time-reversal symmetry the transfer parame¬ 
ters contributing to Eq. (5) satisfy the relations |£ m?n | = 
|£_ m? _ n |. Using these relations, the total exchange pa¬ 
rameter J = J s . s + J s _d is obtained as 


J = ~ 


(2 2 
\ih~2 + u^i 


7 mEs i nEs2 


2 


(6) 


Despite its similar form to in Eq. (4), the above 


TABLE I. The exchange coupling parameters J , and 

Js-d (cm -1 ) for strongly axial magnetic complexes (Fig. 2). J 
corresponds to experimentally extracted exchange parameter. 


System 

Ln 

Ref. 

J 

Js-s 

Js-d 

(a) 

Tb 

[25] a 

b- 

iO 

CO 

1 

-3.58 

0.01 

(a) 

Dy 

[25] B 

-2.97 

-2.51 

-0.46 

(a) 

Ho 

[25] a 

-3.22 

-1.84 

-1.38 

(b) 

Dy 

[26] 

-2.20 

-1.78 

-0.42 

( c ) 

Dy 

[27] 

-0.51 

-0.51 

0.00 


a The exchange parameter J was obtained from experimental 
Ising parameter after extracting the magnetic dipole 
interaction. 


expression involves different orbitals in the second sum¬ 
mation. 

The kinetic exchange interaction between a strongly 
axial doublet and an isotropic spin is described by the 
same Ising Hamiltonian (3) in which one of the pseu¬ 
dospin operators is replaced by the real spin projection 
S z of the corresponding site. The expressions for the 
exchange parameters coincide with Eqs. (4) and (5), in 
which the second summation runs over real orbitals (a) 
for isotropic spin site. Applying the same argument as 
for the orbital m = 0 in the previous case [24], we come 
to the relations |£_ n?a | = |£ n?a | which cancel the terms in 
each bracket of Eq. (5). Thus no s-d kinetic mechanism 
is expected in this case. 

Assessment of J s . s and J s .d in lanthanide complexes 
—. To assess the importance of s-d contribution to the 
exchange interaction in real complexes, we performed a 
density functional theory (DFT) based analysis of J s . s 
and Js_ d for several previously investigated lanthanide 
complexes [25-27]. To this end, we first made the lo¬ 
calization of Kohn-Sham orbitals on the metal centers 
and the bridging ligand. This allowed us to extract the 
metal-ligand transfer parameters (the metal-metal ones 
turned out to be negligibly small in this approach). This 
tight-binding model together with the Hubbard repulsion 
energy (described by one single parameter V due to the 
equivalence of the metal sites, see Fig. 2) was downfolded 
on the ground spin-orbit doublet states of Ln ions (Fig. 
la). This allowed us to calculate straightforwardly the 
energies of ferromagnetic and antiferromagnetic configu¬ 
rations (Fig. lb) and to obtain the corresponding total 
J in Eq. (3). Then, repeating this procedure by blocking 
electron transfer processes from double occupied orbitals 
on the Ln sites, we obtain a net s-s contribution to the 
exchange coupling, J7 S _ S , and finally the s-d contribution: 
ffs-d — J — Js-s- I n this calculations, the parameter U 
was chosen to reproduce the experimental exchange pa¬ 
rameter J (Table I). 

The obtained s-s and s-d contributions are given in 
Table I. The Dy and Ho complexes from isostructural 
series (a) and the Dy complex (b) show that the s-d con¬ 
tribution is by far not negligible in comparison with the 
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(a) (b) (c) 



FIG. 2. The structures of investigated binuclear lanthanide complexes Ln 2 . (a) Ln = Tb, Dy, Ho [25], (b) Ln = Dy [26], and 
(c) Ln =s Dy [27]. Color legend: Ln purple, O red, C gray, N blue, S yellow, Cr orange, Si beige. The pink dashed line is the 
direction of the main magnetic axis and the green arrow is the magnetic moment on Ln ions calculated ab initio [25-27]. 


s-s contribution. The increase of the s-d contribution 
with Ln atomic numbers in the isostructural series (a) is 
explained by the increase of the number of the double- 
occupied orbitals. On the other hand, the vanishing s-d 
contribution in the complex (c) is due to the cancellation 
of the ferromagnetic and antiferromagnetic contributions 
in the expression for J s _a, Eq. (5). 


C. Exchange interaction for non-collinear doublets 

In non-collinear magnetic systems the main magnetic 
axes on sites make an angle </> (Fig. 3a). The new feature 
which appears in this case is that electron can transfer 
to an orbital of a neighbor site in both ferro and anti- 
ferro configurations (see the definition in Fig. 3), with 
the probability depending on <fi. As in the collinear case, 
the single-electron transfer processes cannot switch the 
multiplet components, |J, J) <=* |J, — J), when J is suf¬ 
ficiently large [22]. Therefore, the exchange interaction 
will be described by the same Ising Hamiltonian (3) with 
the difference that now pseudospin operators describe 
momentum projections along corresponding main mag¬ 
netic axes [z\ and z 2 in Fig. 3a). The exchange param¬ 
eter corresponding to s-s processes is obtained as 

= ' Gt + ft) 

X F F ( COs2 |l*m,-n| 2 -sin 2 ^m,„| 2 ) ,(7) 

mE-si n£s2 ' ' 

and contains now both ferro and antiferro contributions. 
On the other hand the exchange parameter for the s- 
d processes remains unchanged, Eq. (5). One can see 
from Eq. (7) that J s . s is not proportional to cos (j) unless 
we have an additional condition |£ m?n | = As 

was discussed above, the latter is fulfilled for interacting 
axial doublet and isotropic spin, in which case also the 


s-d contribution, Eq. (5), vanishes. One should note that 
the transfer parameters in Eqs. (5) and (7) are defined 
for orbitals quantized along main magnetic axes on the 
corresponding metal sites and, therefore, are implicitly 
dependent on angle (j). 


D. Ferromagnetic kinetic exchange interaction 

Contrary to collinear case for which J is always anti¬ 
ferromagnetic (Eq. ( 6 )), the contributions and 
can be of either sign in non-collinear systems, so that 
the resulting exchange interaction can be both ferro and 
antiferromagnetic. To investigate this situation we con¬ 
sider a symmetric homonuclear dimer model. We assume 
that the electron transfer only takes place between one 
pair of orbitals, t = 7 ^ 0 , where ±/x are or¬ 

bital momentum projections on the common axis 2 con¬ 
necting the metals. Figure 3(a) shows calculated J for 
li — 3 as function of </>. For small angles, 0 < </> < 7 r/ 4 , 
J is always antiferromagnetic. In this domain \J\ de¬ 
creases with increasing cf) and at some critical (j) c < 7r/2 
becomes ferromagnetic. Remarkably, the magnitude of 
the ferromagnetic J is of the order ~ t 2 /U and its rel¬ 
ative strength gradually increases when approaching the 
end of the lanthanide series. Figures 3(b),(c) show the 
evolution of J s . s (7) and J s _& (5). For </> < 7 r/ 2 , the s-s 
and s-d processes stabilize the antiferro and ferromag¬ 
netic states, respectively. They are found in competition 
and the latter (s-d) begins to exceed the former (s-s) at a 
critical </> c , which has a simple explanation. The number 
of single-occupied orbitals decreases with the increase of 
the number of / electrons (TV). Following this trend, the 

s will always decrease with N. On the contrary, 
roughly depends on the multiplication of the number of 
single- and double-occupied orbitals. This is the reason 
why it first increases with N till N = 11 and then be¬ 
gins to decrease. As a result, the critical <p c decreases 
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FIG. 3. s-s and s-d exchange contributions in non-collinear system, (a) Definition of ferro (i) and antiferro (ii) ordering for 
non-collinear case. The main magnetic axes z\ and Z 2 are generally non-coplanar. (b) J, (c) J7s- s , and (d) J s .a for Ln 3+ dimers, 
in the symmetric exchange model (inset of plot (a)) as function of the angle cf between the local main magnetic axes. The only 
non-zero transfer parameter is t = 


with the increase of N and in the cases of Ho, Er and 
Tm complexes, the critical (j) c is as small as ca 7t/4. The 
reasons given above explain also the decrease of J in 
the domain 0 < 0 < 7t/2 when moving towards the end 
of the lanthanide series (Fig. 3(a)). For f> > 7t/2, the 
s-s and s-d processes tend to stabilize the ferro and an¬ 
tiferromagnetic states, respectively. In this domain the 
contribution from the s-d processes is dominant because 
the contribution from s-s processes gradually decreases 
with (p and becomes completely quenched at 0 = tt. 

The change of the sign of kinetic exchange parameter 
is not specific only to lanthanides. Similar results are ob¬ 
tained for compounds with transition metal sites in axial 
ground doublet states with unquenched orbital momen¬ 
tum. We obtain again that in the domain 7r/4 < <j) < tt/2 
the exchange parameter becomes ferromagnetic. More¬ 
over, for d 7 and d 8 metal ions this can attain values of 
~ t 2 /U, which corresponds to a very strong ferromag¬ 
netic coupling for transition metal compounds (see Sup¬ 
plemental Materials). 

Effect of non-collinearity on exchange coupling in Er 2 
complex —. The evolution of the exchange parameter 
in function of the angle between the magnetic axes on 
metal sites is studied on the example of an Er 2 com¬ 
plex [28] (Fig. 4a). The calculation has been done in 


full analogy with the previous case of collinear Ln 2 com¬ 
plexes (Fig. 2 and Table I). Similarly to the model cal¬ 
culations (Fig. 3), with the increase of </> the antifer¬ 
romagnetic exchange interaction becomes ferromagnetic 
around (f> c « 27r/5 (Fig. 4b). The shift of (f c in compar¬ 
ison with the model calculations is due to the existence 
of many electron transfer processes in this complex. In 
real systems, the direction of main magnetic axes could 
be controlled by varying ligand environment. 


III. CONCLUSIONS 

In this work, we investigated the kinetic exchange 
interaction between axial magnetic doublets with un¬ 
quenched orbital momentum. We find a new mecha¬ 
nism of exchange interaction based on electron transfer 
between single- and double-occupied orbitals. Contrary 
to conventional spin systems, the s-d kinetic contribu¬ 
tion found here is not related to Goodenough’s mecha¬ 
nism (2), arising due to the Hund’s rule coupling (Jh) 
on metal sites, but due to the second-order kinetic mech¬ 
anism (1). On this reason, this kinetic contribution is 
as strong as the conventional kinetic exchange between 
single-occupied orbitals but, at variance with the lat- 
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FIG. 4. Ferromagnetic exchange interaction induced by strong magnetic anisotropy, (a) The definition of the angles cj) and if; 
defining the relative orientation in the Er 2 complex with equivalent metal sites, (b) Variation of J with respect to 0 and i/;. 
The internal dashed line corresponds to the angle 7r/4 and the external to the angle 7r/2 between the magnetic axes. The blue 
and the red regions stand for the negative and positive values of J, respectively (cm -1 ). 


ter, can be ferromagnetic. In non-collinear magnetic 
systems the s-d kinetic mechanism can cause an over¬ 
all ferromagnetic exchange interaction of the order of 
starting from angles ~ 7t/4 between main magnetic 
axes. These conclusions are fully supported by quan¬ 
tum chemistry based analysis of Ln 2 complexes. The 
key feature underlying the new mechanism is that the 
double-occupied orbitals change under time inversion in 
strongly anisotropic sites due to unquenched orbital mo¬ 
mentum. This is found in sharp contrast to the case of 
isotropic and weakly anisotropic sites, where no change 
of double-occupied orbitals occur under time inversion. 
The obtained results offer a new view on the exchange 
interaction in lanthanides, actinides and transition metal 
ions with unquenched orbital momentum. In particular, 
they show the way to achieve strong ferromagnetic cou¬ 
pling between metal ions, a long sought goal in magnetic 
materials [4]. 

MATERIALS AND METHODS 

The DFT calculations have been done with the ORCA 
package [29], using B3LYP exchange-correlation func¬ 


tional [30], in which the Hartree-Fock contribution to 
the exchange part was increased from 20 % to 40 %. 
This was done to reproduce the experimental isotropic 
exchange parameters in isostructural Gd 2 analogues of 
investigated complexes. The derivation of tight-binding 
Hamiltonian for localized Kohn-Sham orbitals and the 
projection of the Hubbard model on the ground doublets 
of investigated Ln 2 complexes is described in Supplemen¬ 
tal Materials. 
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“New mechanism of kinetic exchange interaction 
induced by strong magnetic anisotropy” 

This material contains: 

1) The lanthanide complexes and modeling for the elec¬ 
tronic structures; 

2) The density functional theory calculations; 

3) The extraction of the Ising exchange parameters; 

4) The s-s and s-d exchange mechanisms in the d metal 
complexes. 


The model Hamiltonian Eq. (SI) is not the same as 
the model Hamiltonian used for the derivation of the 
exchange interactions in the main text. In the latter 
model Hamiltonian, the electron transfer between the 
metal sites and ligand is reduced to the transfer parame¬ 
ter between the metals. However, in the DFT based cal¬ 
culations of the exchange parameters, such simplification 
is not always possible because the / orbital levels and 
ligand levels are sometimes close to each other. Thus, 
we take Eq. (S3) for our calculations of the lanthanide 
complexes. 


I. MODEL HAMILTONIAN FOR BINUCLEAR 
MAGNETIC SYSTEMS 

As the anisotropic magnetic systems, we treat binu- 
clear lanthanide complexes [S1-S4]. In order to describe 
the low-energy exchange states of these systems, the Hub¬ 
bard Hamiltonian was used: 

H = H 0 + H t + H bu (SI) 

Hq — ^ ^ e r n ra ', (S2) 

rcr' 

Ht = E EE (r£r>L(^)4n*<w 

i= 1,2 mcr rcr' 

+ r r L ;d;(fii)4 r ^w), (S3) 

y ' y ; ■ (S4) 

i= 1,2 mcr^m'cr' 

Here, Ho, H t , H b [ indicate the orbital energy level of the 
bridging atoms with respect to the / orbital level, the 
electron transfer Hamiltonian between metal site and the 
bridging atoms (ligand), and H b [ is the Coulomb repul¬ 
sion on / metal sites, respectively. In the Hamiltonian, i 
(= 1,2) is the metal (lanthanide) center, L is the bridging 
ligand, m (= 3, 2,..., —3) is the projection of the orbital 
angular momentum of / atomic orbital, r is the orbital 
energy level of the ligand, a, a' is the spin projection, 
Cima (cJmcr) an d ra (4,ra) are t ^ ie annihilation (cre¬ 
ation) operators on the ma orbital of metal site and ra 
orbital of the ligand L, respectively, h irn(7 = cJ ma c im(T 
and fiL r( j = c\ jrcr CL rcr are the number operators, r^ r 
(= t^*) is the transfer parameter between magnetic cen¬ 
ter and ligand, and Uf is the Coulomb repulsion energy 
on / metal site. Since the transfer parameter between 
metal sites are small, we neglect the direct electron trans¬ 
fer between the metal centers. The Coulomb repulsion 
on the ligand and the intersite one between the metal 
center and the ligand as well as the Hund’s rule cou¬ 
pling are neglected because they are smaller than Uf. 
Based on the Hamiltonian, we describe the low-energy 
magnetic states. The transfer parameters are extracted 
from density functional theory (DFT) calculations and 
the Coulomb repulsions are determined to reproduce the 
experimental exchange parameters. 


II. DFT CALCULATIONS 

In order to find a density functional which is suit¬ 
able for the electronic structure calculations of lanthanide 
complexes, we calculated (i) the exchange interactions 
and (ii) the binding energies with several methods (DFT 
and Hartree-Fock (HF) method), and compared the re¬ 
sults with experimental data. We found that the hybrid 
functional including about 40 % of HF exchange is suit¬ 
able to express one-orbital parameters (e r and r^ r in Eqs. 
(S2) and (S3)) of lanthanide complexes. For the quan¬ 
tum chemistry calculations, we used ORCA3.0 package 
[S5]. 


A. Isotropic exchange parameters for Gd 3+ 
complexes 


In general, the spin-orbit coupled ground states of the 
lanthanide ion (Dy 3+ , Er 3+ , etc.) cannot be adequately 
treated within single Slater determinant approach. On 
the other hand, the ground state of half-filled Gd 3+ ion 
is in a good approximation described by pure spin state, 
and the state with maximal spin projection is a single 
Slater determinant. Thus, for the quantum chemistry 
calculations, we replaced the lanthanide ions (Dy 3+ and 
Er 3 + ions) by Gd 3+ ions. As the functional, we chose 
B3LYP with various HF exchange contributions (20 - 
60 %). The variation of the contribution of the HF 
exchange is to reduce the exaggerated electron transfer 
and to increase the underestimated Coulomb repulsion 
on site. The molecular structures were taken from the 
x-ray diffraction data. As the basis set, SVP was used. 

The exchange parameter between metal sites was ob¬ 
tained within broken-symmetry approach [S6]. The ex¬ 
change parameter J defined by Heisenberg Hamiltonian, 

H = - JSi • S 2 , (S5) 


is estimated using the high-spin and broken-symmetry 
(low-spin) states: 


j = ^ (~ Ehs ~ ~ Ebs ^ 

(S 2 )hs — (S 2 )bs 


(S6) 


Here, Fhs and E^s are the ground electronic energies 
obtained from the high-spin (HS) and broken-symmetry 
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(BS) calculations, and (S 2 )hs and (S 2 )bs are the expec¬ 
tation values of the magnitude of the total spin. 

The obtained isotropic exchange parameters and those 
from experimental data are tabulated in Table SI. We 
find B3LYP functionals with 40-50 % of Hartree-Fock 
exchange well reproduce the the experimental ones. 


B. Binding energy vs. photoemission 
measurements 

The binding energies of tris-cyclopentadienyl lutetium 
(LuCps) were estimated and compared with experimen¬ 
tal data. The experimental binding energies are taken 
from the photoelectron spectra (PES) of LuCp 3 in gas 
phase [S8]. We chose Lu 3+ ion complex because the 4/ 
orbitals are completely filled in the ground state, which 
can be described by singe Slater determinant. For the 
calculations, we used the x-ray diffraction structure of 
CeCp 3 [S9] as the structure of LuCp 3 since their struc¬ 
tures are similar to each other and the energy scale of the 
PES is much larger than the change in energy due to the 
structure. The binding energies were estimated applying 
the Koopmans’ theorem. 

The calculated binding energies are shown in Table S2. 
In comparison with the experimental data, the B3LYP 
(20 %) and the HF calculations underestimates and over¬ 
estimates the 4/ binding energies, respectively. Better 
agreement is obtained when the contribution of the HF 
exchange is about 40 %. Since both the exchange inter¬ 
action and the binding energies obtained with the hybrid 
B3LYP functional with about 40 % of the HF exchange 
are close to the experimental data, we used the functional 
for the calculations of the transfer Hamiltonian. 


III. DERIVATION OF ISING EXCHANGE 
PARAMETERS FOR LN 3+ COMPLEXES 

In order to project the electronic states from the DFT 
calculations into the model Hamiltonian (SI), we local¬ 
ized the valence Kohn-Sham orbitals (Pipek-Mezey lo¬ 
calization [S10]). The Kohn-Sham orbitals are divided 
them into the magnetic 4/ orbitals, the bridging ligand 
orbitals, and the other ligand orbitals, where the last 
orbitals are not important for the exchange interaction 
from each other. The localization on the magnetic core 
was measured by the Mulliken population on the / type 
orbitals of the metal centers and the p type orbitals of the 
bridging ligand. With the use of the localized orbitals as 
the basis set, we expressed the Kohn-Sham Hamiltonian 
matrix hxs for the magnetic core. Here, the localized 
/ orbitals were transformed into the eigenstates of the 
atomic orbital angular momentum l z whose quantization 
axes agree with the ab initio main magnetic axes [S1-S3], 

| ima). On the other hand, the ligand orbitals | Lra) are 
chosen so that the bridging ligand part of Iiks becomes 
diagonal. Then, the diagonal elements of the Hamilto¬ 


nian matrix hxs are used as the orbital energy levels 
(the / orbital level is averaged) in Eq. (S2), and the 
elements in the off-diagonal block matrices between the 
metals and the bridging ligand as the transfer parameters 
in Eq. (S3). 

The matrix of model Hamiltonian (SI) was calculated 
using the ground ferromagnetic, the ground antiferro¬ 
magnetic electron configurations (Fig. lb in the main 
text), and the configurations with one-electron transfer 
between metal sites and between metal and ligand as 
the basis set. Diagonalizing the Hamiltonian matrix, we 
obtain the ferromagnetic and antiferromagnetic states. 
When there is no energy level between these two states, 
we can project these states into Ising Hamiltonian (Eq. 
(3) in the main text). The Coulomb repulsion Uf is deter¬ 
mined to reproduce the experimental exchange parame¬ 
ter. 

The obtained exchange parameters are shown in Ta¬ 
ble I in the main text. The Coulomb repulsion energies 
are 1.50, 1.40, 1.20 eV for the Tb, Dy, Ho complexes in 
the series (a), 1.48 and 1.44 eV for the complexes (b) 
and (c), respectively. In all cases, there is no level be¬ 
tween the ground ferro- and antiferromagnetic levels. At 
a glance, Uf looks too small (cf. Ref. SI 1), whereas it 
is not too small because the effect of the Coulomb re¬ 
pulsion is partly included in the Kohn-Sham orbital en¬ 
ergies within the mean-field approximation. Thus, Uf 
is the difference between the Hubbard U and the mean- 
field value. The ferromagnetic and the antiferromagnetic 
ground states are mainly contributed by the electron con¬ 
figurations without electron transfer (the type of Fig. lb 
in the main text). Their contributions (probabilities) to 
the ground state are 98.3 %, 99.4 %, 98.9 % for the Tb, 
Dy, Ho complexes of the series (a), 99.3 % for the com¬ 
plex (b), 99.7 % for the complex (c), and 97 % for the 
Er complex. These high probabilities guarantee the va¬ 
lidity of the description of the low-lying states by the 
pseudospin Hamiltonian. 


IV. S-S AND S-D EXCHANGE MECHANISMS 
IN d METAL COMPLEXES 

Within the simple two-sites model for the non-collinear 
doublets in the main text, we calculated the exchange 
interaction parameter for d metal complexes assuming 
that the ground state is | J, d= J). The exchange parameter 
J for the d 6 , d 7 and d 8 ions are shown in Fig. SI (a). As 
in the case of the / metal ions, the exchange becomes 
ferromagnetic as the increase of angle <fr. 

Although the crystal field in transition metal complex 
is stronger than in lanthanide complex, the crystal field 
level with unquenched orbital can be obtained with suit¬ 
able symmetry of the ligands. For example, axial (d 6 , d 7 ) 
or trigonal (d 8 ) crystal field splits the d levels into two 
doublets and one nondegenerate state [S12]. With the 
splitting, the ground doublets of the d 6 metal originates 
from the J-multiplet due to the strong Hund’s rule cou- 
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TABLE SI. The exchange parameters extracted from the broken-symmetry DFT calculations and experimental magnetic 
susceptibility (cm -1 ). The numbers in the first row indicate the contribution of the HF part in the exchange-correlation 
functional. The last column are the experimental data (Exp.) extracted using the Lines model or Heisenberg model. The 
experimental data of the Gd complexes (a), (b), (d) are taken from Refs. SI, S2, and S7, respectively. The effect of the 
magnetic dipolar interaction is removed only for the complex (a). 



20 % 

25 % 

30 % 

35 % 

40 % 

45 % 

50 % 

55 % 

60 % 

Exp. 

(a) 

-0.300 

-0.268 

-0.235 

-0.213 

-0.194 

-0.180 

-0.160 

-0.150 

-0.133 

-0.17 [SI] 

(b) 

-0.319 

-0.280 

-0.249 

-0.231 

-0.207 

-0.190 

-0.173 

-0.157 

-0.142 

-0.210 [S2] 

(d) 

-1.468 

-1.057 

-0.785 

-0.603 

-0.478 

-0.385 

-0.312 

-0.260 

-0.222 

-0.448 [S7] 


TABLE S2. Binding energies (eV). 4/, Cp (cr), and Cp (n) are the nature of the Kohn-Sham orbitals for the corresponding 
range of energies. The experimental energies indicate the peaks of the photoelectron spectra [S8]. 




B3LYP 




20 % 

25 % 

30 % 

35 % 

40 % 

45 % 

4/ 

12.47 - 12.59 

13.21 - 13.33 

13.96 - 14.08 

14.71 - 14.82 

15.46 - 15.58 

16.21 - 16.33 

Cp (cr) 

9.97 - 11.22 

10.26 - 11.59 

10.56 - 11.95 

10.85 - 12.32 

11.14 - 12.68 

11.43 - 13.05 

Cp (tt) 

5.24 - 6.71 

5.42 - 6.92 

5.60 - 7.12 

5.78 - 7.33 

5.97 - 7.54 

6.15 - 7.74 



B3LYP 


HF 

Exp. 



50 % 

55 % 

60 % 




4/ 

16.97 - 17.09 

18.04 - 17.68 

18.63 - 18.54 

23.29 - 23.21 

14.30, 15.74 


Cp (cr) 

11.73 - 13.42 

12.02 - 13.79 

12.31 - 18.35 

15.92 - 13.77 

12.4 


Cp (tt) 

6.33 - 7.95 

6.52 - 8.16 

6.70 - 8.37 

9.13 - 7.22 

7.28, 8.81 



pling (Fig. Sl(b)). In the case of the d 7 and d 8 systems, 
the orbital momentum is unquenched in the presence of 


the spin-orbit coupling which exceeds the crystal field 
splitting of the orbital levels (A in Fig. Sl(b)), which 
could be observed even in 3d metals systems. 
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FIG. SI. (a) J for several d-metal complexes as function of the angle between local main magnetic axes, (b) d 6 , d 7 and 
d 8 electron configurations corresponding to axial doublets with maximal Mj. The non-aufbau population in the last two 
electron configurations is achieved for relatively small separation A between orbital levels with non-zero m (compared to spin- 
orbit coupling). This order of orbitals is realized, for example, in crystal fields with axial symmetry (d 6 , d 7 ) and/or trigonal 
symmetry (d 8 ) [SI2]. 
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